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Abstract 

We consider an extension of the Standard Model by a singlet scalar 
that accounts for the dark matter of the Universe. Within this model 
we compute the expected gamma ray flux from the annihilation of dark 
matter particles in a consistent way. To do so, an updated analysis of 
the parameter space of the model is first presented. By enforcing the 
relic density constraint from the very beginning, the viable parameter 
space gets reduced to just two variables: the singlet mass and the 
higgs mass. Current direct detection constraints are then found to 
require a singlet mass larger than 50 GeV. Finally, we compute the 
gamma ray flux and annihilation cross section and show that a large 
fraction of the viable parameter space lies within the sensitivity of 
Fermi-GLAST. 



1 Motivation 



A simple extension of the Standard Model that can explain the dark matter 
is the addition of a real scalar singlet and an unbroken Z2 symmetry under 
which the singlet is odd while all other fields are even. Such a singlet, which 
couples directly only to the higgs boson and to itself, may indeed have the 
right relic density -in the standard cosmological model- to explain the ob- 
served dark matter abundance. Even though this singlet extension has been 
studied several times [H [21 El IH E] , a complete and consistent computation 
of the expected gamma ray flux from the annihilation of singlet dark matter 
has yet to be published. In this paper we will fill that void. First, an updated 
analysis of the parameter space of the model will be presented. Then, we 
compute the direct detection cross section and show that present constraints 
require a singlet mass larger than 50 GeV. Finally, after obtaining the gamma 
ray flux and total annihilation cross section, we show that over most of the 
viable parameter space the singlet scalar model of dark matter is detectable 
by Fermi-GLAST [B]. 

The observation of gamma rays originating in the annihilation of dark 
matter particles is one of the most promising avenues to determine the nature 
of dark matter. The recently launched Fermi Gamma-ray Space Telescope 
(FGST), formerly GLAST will improve, with respect to its predecessor 
EGRET [7], the sensitivity to gamma rays in the energy range between 20 
MeV and 10 GeV by more than one order of magnitude. More importantly, 
FGST will extend the high energy range to about 300 GeV, making it an ideal 
experiment to search for gamma rays from WIMP dark matter annihilation. 
It is therefore critical to determine the expected gamma ray flux within 
diverse scenarios accounting for the dark matter. 

The lightest neutralino in supersymmetric models is by far the most com- 
mon dark matter candidate examined in the literature. Even if less com- 
pelling from a theoretical point of view, the singlet scalar considered in this 
paper offers an interesting alternative to neutralino dark matter. Depending 
only on two new parameters, the singlet scalar model is very predictive and 
could be easily falsified. Moreover, the dark matter candidate is a scalar 
instead of a Majorana fermion and its main annihilation channels do not co- 
incide with those of the neutralino. Direct and indirect detection signals are 
thus expected to be different. Finally, the disparity between the implications 
for colliders searches of these two dark matter models could not be more 
marked. The singlet scalar model predicts the existence of one additional 
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degree of freedom rather than a full spectrum of superpartners. The LHC, 
therefore, may soon shed some light on the identity of the dark matter parti- 
cle and, in particular, on its supersymmetric or non-supersymmetric nature. 
In the meantime, it is important not to restrict ourselves to supersymmetric 
candidates. 

Over the years, several authors have studied the phenomenology of the 
singlet extension of the Standard Model. The singlet scalar as a dark mat- 
ter candidate was initially proposed by McDonald [Ij and was subsequently 
analyzed by Burgess et. al [2J. Recently, Barger et al. [1] considered models 
with and without the Z2 symmetry and investigated their expected phe- 
nomenology at the LHC. Regarding dark matter, they computed the singlet 
relic density as well as its direct detection cross section. Later on, in [5], 
the gamma ray flux was computed but only for two specific values of the 
singlet-higgs coupling. In this paper we update and expand these previous 
results in several respects. First, we consider the full mass range, including 
nis < M]Y and ms > rrih, for the singlet scalar. Second, we use the precise 
determination of the dark matter density obtained by the WMAP experiment 
[S] as well as the accurate computation of the relic density by micrOMEGAs 
[9] -which includes all tree-level annihilation processes- to obtain the viable 
parameter space of the model. Third, we compute the direct detection cross 
section and take into account the constraints from current experiments. Fi- 
nally, we use state of the art techniques, as implemented in micrOMEGAs 
[Hj, to compute the expected gamma ray flux along the viable regions of the 
singlet model. 

The paper is organized as follows. In the next section we will introduce 
the model Lagrangian and will identify the new parameters that it contains. 
In section [3] the relic density is computed and used to obtain the viable 
parameter space. Direct detection rates are then calculated and compared 
with the sensitivity of present and planned experiments in section HI Finally, 
in section [5l we compute the expected gamma ray flux and annihilation cross 
section and show that Fermi-GLAST will probe a large fraction of the viable 
parameter space of the singlet scalar model. 
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2 The model 



The Lagrangian that describes the model with an additional scalar singlet, 
S, is 

L = CsM + Id.Sd'^S - - ^S' - XS'H^H , (1) 

where Csm denotes the Standard Model Lagrangian and H is the higgs dou- 
blet. This Lagrangian is the most general renormalizable one that is compat- 
ible with the SU{3) x SU{2) x U{1) gauge invariance and with the symmetry 
S —>■ —S. The scalar singlet extension of the standard model, therefore, 
contains only 3 new parameters: mg. A, and A5. Because it only determines 
the strength of the singlet self-interactions, A5 is unconstrained and largely 
irrelevant to the phenomenology of the model. In the following we will simply 
require A5 < 1 so as to guarantee a perturbative treatment. Notice, from 
dH), that the singlet couples to Standard Model fields only through the higgs 
boson and that such interaction is determined by the parameter A. 

A detailed analysis of the scalar potential was already presented in [2]. 
Here, we briefly review, for completeness, the constraints that the potential 
must satisfy. In the unitary gauge, the scalar potential takes the form 

V = ^S^ + \s^h^ + ^S^ + ^{h^ - vl^f (2) 

where vew = 246 GeV and \h is the higgs quartic coupling. The configu- 
ration = and 7^ is a local minimum of V provided that v\y^ > 
and ml + Xv%^r > 0. Another local minimum, with h = and = —ml/X, 
will exist whenever ml < and —Xm^ > XsXhV%^r. In such case, to en- 
sure that the former is the potential's global minimum we must require that 
< —ml < v]^yyy/XhXs- Once these conditions are satisfied, the S-dependent 
part of the potential can be rewritten as 

V = im|^2 + ^S' + XvEwS'h + ^S'e (3) 

where m| = ml + Xvj^y^r, and h represents the physical higgs boson with mass 
m\ = XhV^^^. In the following we take ms and A as the free parameters of 
the singlet scalar model. 
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Figure 1: The dark matter density as a function of ms for A = 0.1 and 
different values of the higgs mass. 



3 The viable parameter space 

In this section we compute the rehc density of the singlet scalar and use it 
to impose the dark matter constraint, Vlsh^ = 0.11 [8]. From it, we obtain 
the viable parameter space of the model with a singlet scalar. 

Singlets can annihilate through s-channel higgs boson exchange into a 
variety of final states: //, W^W~, Z^Z^, and hh. Additionally, they can 
also annihilate into hh either directly or through singlet exchange. As a 
general rule, the final state W^W~ tends to dominate the total annihilation 
cross section whenever such channel is open. A light singlet, ms < My/, 
will annihilate mainly into the hh final state. An intermediate mass singlet, 
rriw < ms < mt, annihilates mostly into W~^W~, with additional contribu- 
tions from Z^Z^ and, if allowed, hh. For a heavier singlet, ms > mt, the 
pattern is similar, as the tt channel gives a non-negligible but subdominant 
contribution. To accurately compute flsh"^ we use the micrOMEGAs package 
[9], which can calculate the relic density in a generic dark matter model. 

Figure [1] shows the relic density as a function of ms for A = 0.1 and differ- 
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Figure 2: The viable parameter space of the scalar singlet model. Along the 
lines the dark matter constrained is satisfied. We use mh = 120, 150, 180, 200 
GeV as reference values for the higgs mass. The grey area surrounding the 
line corresponding to m/j = 120 GeV shows the region compatible with the 
observed dark matter density at 2a. 



ent values of the higgs mass. Notice that the scalar singlet model can explain 
the dark matter naturally -that is, without any fine-tuning in the parame- 
ters. Indeed, for A = 0.1 and nis around the electroweak scale, the predicted 
relic density lies in the correct range to be compatible with the observations. 
The most noticeable feature from this figure is the drastic suppression of 
the relic density that takes place at the higgs resonance. In fact, around 
2ms ~ "^h the relic density is orders of magnitude smaller than anywhere 
else. The effect of the threshold is also seen to be important. Above it, 
ms > mw, the W~^W~ annihilation channel is open and consequently the 
relic density tends to be smaller than below it. 

It is clearly seen from figure [1] that the value of the higgs mass, even 
if a standard model parameter, is critical for the computation of the sin- 
glet relic density. From direct searches at LEP, a lower limit on the higgs 
mass can be obtained, rrih > 114.4 GeV [lOj. When this bound is com- 
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bined with electroweak precision measurements, an upper limit -at the 95% 
C. L.- of 182 GeV [TT]| is derived. Throughout this paper we will use 
rrih = 120, 150, 180, 200 GeV as reference values for the higgs mass. 

For any given pair {mh,ms) there exists a unique value of A such that 
the dark matter constraint, Qsh"^ = 0.11, is fulfilled. By imposing the dark 
matter constraint, therefore, the variable A can be effectively eliminated for 
given values of rrih and ms, reducing the viable parameter space to a two 
dimensional volume. In figure [2] we show, in the plane {X,ms), lines that are 
compatible with the observed dark matter density for different values of ruh- 
Away from the higgs resonance, the typical value of A is O{10~^ — 10~^). 
Light singlets, which annihilate mostly into bb, require larger values of A 
to obtain the observed relic density. In contrast, at the higgs resonance, 
the annihilation tends to be more efficient and much smaller values of A are 
needed to satisfy the dark matter constraint. To be concrete, we will not 
consider singlet masses above 300 GeV, though such high values are also 
allowed. Figure [2] defines what we call the viable parameter space of the 
singlet scalar model. In the following we will compute the direct detection 
rates and the gamma ray flux only along these viable lines. 

4 Direct detection 

Dark matter particles can scatter elastically on nuclei and be detected through 
nuclear recoil in direct detection experiments. In the non-relativistic limit, 
the dark matter-nucleon amplitude receives two contributions: the scalar 
or spin- independent interaction and the axial- vector or spin-dependent part. 
For the singlet scalar, the spin- dependent interaction vanishes so only the 
spin-independent part can give a signal. 

Figure [3] shows the spin-independent proton-singlet cross section as a 
function of the singlet scalar mass (see [9] to find out how this cross sec- 
tion is computed). Since the higgs- singlets coupling is small close to the 
higgs resonance, see figure [21 the cross section is highly suppressed in that 
region. A heavy singlet, ms ^ 150 GeV, has an interaction cross section 
around 10^^ pb. Because the singlet interacts with nucleons via t-channel 
higgs exchange, the cross section typically decreases with the higgs mass, 
as observed in the figure. For reference, current constraints from XenonlO 
[T2] and CDMS [13j are also displayed. They rule out singlet masses below 
50 GeV independently of the higgs mass. Nevertheless, the possibility of ex- 
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Figure 3: The spin-independent proton-singlet cross section as a function of 
rris for different values of the higgs mass. The thin lines show the present 
constraint from XENONIO and CDMS. The dotted line corresponds to the 
expected sensitivity of SuperCDMS. Along the lines flsh"^ = 0.11. 



plaining the DAMA signal [H] with a light singlet, ms ~ 5 GeV, was recently 
explored in [15] . Finally, notice from the figure that future experiments, such 
as SuperCDMS [16], will probe a significant region of the viable parameter 
space. 

5 The gamma ray flux 

The role of the indirect detection of dark matter -that is, the detection of dark 
matter annihilation products- is complementary to that of direct detection 
searches and will be crucial in future dark matter studies. In principle, dark 
matter annihilations could be observed through gamma rays, neutrinos, or 
antimatter. Among them, the simplest and more robust is the gamma ray 
signal. 

The gamma ray flux above some energy threshold Et^r from a direction 
forming an angle with respect to the galactic center can be expressed as 
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^{Ethr) =0.94 X lO^^^cm^^s"^ 

(4) 

where J{ip) is the dimensionless hne of sight integral around the direction ip, 
and J{AQ) is the average of Jiip) over the spherical region of solid angle AQ, 

Jim = ^l^^ ^(^)^^ (6) 

From this expression we see that the gamma ray flux depends not only on 
particle physics parameters, such as (av) and ms, but also on the unknown 
distribution of the dark matter, parameterized by J{AQ). 

In micrOMEGAs the photon spectrum from dark matter annihilations is 
obtained with a procedure similar to that implemented in DarkSUSY [T7] . 
Such procedure relies on the use of tables of 7 production -obtained from 
PYTHIA for each of the different basic channels: qq, fi^fJ^^, r^r^, 

W^W~ and ZZ. Concerning the distribution of dark matter, in our analysis 
we take a NFW profile [19], which has a p(r) oc behaviour at small r, 
as the canonical example. For comparison we show, in table [1], the value 
of J{AVL) for three different halo profiles -including NFW- and two typical 
values of Af2. 



Profile J{An = IQ-^ sr) JjAQ = 10'^ sr) 

NFW 1.21 ■ 10^ 1.26 • 10^ 

Moore 1.05 ■ 10^ 9.46 • 10^ 

Modified isothermal 3.03 ■ 10^ 3.03 ■ 10^ 



Table 1: Values of {J{0))^q for two different Af2's and for three different 
density profiles. See [20] for details. 
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Figure 4: The integrated photon flux above 1 GeV as a function of the singlet 
scalar mass. The predicted gamma flux is from a AQ = 10^^ sr region around 
the galactic center for a NFW profile. 



Figure m shows the integrated gamma ray flux from the Galactic Center as 
a function of rris for different values of rrih- All these models are compatible 
with the EGRET constraint [20l [2T] . Notice that the flux has a dip not only 
at the higgs resonance but also at the W threshold. This latter feature is 
due to the effect known as annihilation into forbidden channels [22]. The 
l/m| dependence, from is clearly visible in the figure. Lighter singlets 
generically yield a larger 7 ray flux. 

Even though the Galactic Center is expected to be the strongest source 
of 7 rays from DM annihilation, it is not necessarily the best place to observe 
them. To do that, one must be able to distinguish the DM signal from the 
Galactic diffuse background, a task that, given its large backgrounds, might 
be more difficult at the Galactic Center. In fact, a recent study, [23], found 
that, for a NFW profile, the Fermi-GLAST sensitivity to a 7 ray signal from 
the Galactic Halo (excluding the Center) is larger than that from the Galactic 
Center -compare figures 5 and 9 in [23]. In the following, we will use those 
results, in particular figure 9, to determine the Fermi-GLAST sensitivity to 
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Figure 5: The total annihilation cross section at low velocities as a function 
of ms for different values of rrih. The dashed (dotted) line shows the (av) 
required to observe a DM annihilation signal at 3a significance with one 
year of Fermi-GLAST data considering the conventional (optimized) diffuse 
model as background -see figure 9 in [23\ . 



the singlet scalar model of dark matter. 

In figure [5] we show {av), at small v, as a function of the singlet mass and 
different values of rrih. It is clearly seen from the figure that away from the 
higgs resonance and the W threshold (av) is essentially constant and equal 
to the so-called typical annihilation cross section, (av) ~ 3 x 10~^^cm^s~^. 
The dashed line shows the (av) required to make an observation of DM 
annihilation at 3a significance with one year of FGST data considering a 
conventional diffuse model [21] as background. The dotted lines shows the 
analogous quantity but using the optimized diffuse model [25| as background. 
These two lines were obtained by assuming that the region within 10° of the 
GC is excluded from the analysis. Notice from the figure that for the opti- 
mized model most of region with ms < 175 GeV is "detectable by FGST". 
For the conventional diffuse model, the detectable region extends to much 
higher masses. It is fair to say, therefore, that Fermi-GLAST will probe a 
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significant region of tlie parameter space of the singlet scalar model of dark 
matter. 

6 Conclusion 

We have studied in detail the singlet scalar model of dark matter. First, we 
introduced the model and presented an updated analysis of its parameter 
space. The dark matter constraint was found to reduce the viable parameter 
space to just two variables: the singlet mass and the higgs mass. Then, we 
computed the direct detection cross section and found that singlet masses 
below 50 GeV are already ruled out by the recent data from CDMS and 
XenonlO. Finally, we calculated the expected gamma ray flux from the 
Galactic Center as well as the total annihilation cross section and showed 
that a significant region of the parameter space will be probed by Fermi- 
GLAST. 
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